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Two qubit entanglement can be induced by a quantum data bus interacting with them. In this 
paper, with the quantum spin chain in the transverse field as an illustration of quantum data bus, we 
show that such induced entanglement can be enhanced by the quantum phase transition (QPT) of 
the quantum data bus. We consider two external spins simultaneously coupled to a transverse field 
Ising chain. By adiabatically eliminating the degrees of the chain, the effective coupling between 
these two spins are obtained. The matrix elements of the effective Hamiltonian are expressed in 
terms of dynamical structure factor (DSF) of the chain. The DSF is the Fourier transformation of 
the Green function of Ising chain and can be calculated numerically by a method introduced in [O. 
Derzhko, T. Krokhmalskii, Phys. Rev. B 56, 11659 (1997)]. Since all characteristics of QPT are 
embodied in the DSF, the dynamical evolution of the two external spins displays singularity in the 
vicinity of the critical point. 



PACS numbers: 03.65.Ud, 05.50.1q, 05.70.Jk, 75.10.Pq 

I. INTRODUCTION 

Entanglement lies at the heart of quantum mechanics 
and thus can be regarded as a resource for quantum in- 
formation processing. In recent years, many people have 
demonstrated that quantum entanglement can offer an 
intrinsic clarification of quantum criticality of a many 
body system. For example, in Refs. 1, 2] it was proved 
that the derivative of the nearest-neighbor entanglement 
diverged at the critical point. Furthermore, one of the au- 
thors CPS and his collaborators studied the dynamical 
ultrasensitivity of the induced quantum critical system, 
which concerned the excited states as well as the ground 
state d, 0] . Usually, the former has not been discussed 
in the investigation for the above mentioned intrinsic en- 
tanglement in many body systems. 

In this paper, we consider the quantum entanglement 
of the many body problem in the point view of quantum 
information processing, whether or not quantum critical- 
ity of the many body system (as a quantum data bus) can 
enhance the entanglement of the external qubits interact- 
ing with the quantum data bus. Actually, the phase tran- 
sition indeed creates some entanglement. A case in point 
is the superconducting phenomenon 0, Q . As shown in 
the BCS theory 0, Q, a Cooper pair is created when a 
conductor transits from a normal state to a supercon- 
ducting state. In such a pair, two electrons form a cor- 
related entire over hundreds of nanometers, which is a 
long distance in the microscopic world. We remark that 
the phase transition of superconducting happens at finite 
temperature but this paper will focus on the occurrence 
of quantum phase transition (QPT) at zero temperature. 

QPT is of critical importance to the quantum statisti- 
cal physics. Generally speaking, QPT takes place at zero 
temperature. It's the situation where only uncertainty 
principle plays the major role while the fluctuation due 
to the finite temperature does not. As some parameter 
is varied, a qualitative change occurs in the ground state 



of a quantum many-body system due to QPT [lo| . 
At the critical point, long range correlation also develops 
in the ground state. This long range correlation intu- 
itively exhibits greater quantum entanglement between 
two points. Refs.fll, indeed showed the entanglement 
became singularly longer at the critical point. These con- 
siderations directly motivate us to study the problem. 

On the other hand, the interesting thing of this pa- 
per is its meaning for the detection of QPT by coupling 
the quantum critical system to the external detector-two 
qubits. And various efforts have been devoted to this 
field. In Ref.[3], when an external spin underwent a tran- 
sition from a pure state to a mixed state, the decay of 
the Loschmidt echo of its coupling environment described 
by a transverse field Ising model (TFIM) was greatly en- 
hanced. Others suggested two external spins to detect 
QPT by simultaneously coupling to an XY model en- 
vironment [ill , [l^ . The former researches showed the 
critical phenomena of QPT with exact solvable models. 
However, they mainly focused on a^a'' interaction be- 
tween the external spins and the the Ising chain. Thus, 
by extending z-z coupling to a more general coupling 
form, we explore the possibility of detecting QPT by two 
central spins. 

The rest of paper is organized as follows. In the next 
section, we describe the model as two external spins cou- 
pled to a transverse field Ising chain. The calculation of 
the effective Hamiltonian between the two spins is out- 
lined using Frohlich transformation in Sec. HI. The ma- 
trix elements of the effective Hamiltonian are given in 
terms of dynamical structure factor (DSF). In Sec. IV, 
the DSF is numerically calculated and the relation be- 
tween the coupling constants of two spins and the vari- 
able parameter is given. Finally, the significant results 
are concluded in Sec. V. 
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FIG. 1: (color online) Schematic diagram of two external spins 
(C) simultaneously coupling to a 1-D Ising chain (E). 



II. MODEL DESCRIPTION 

Two external spins are simultaneously coupled to an 
environment described by a one dimensional transverse 
field Ising chain. By exchanging spin angular momentum 
with the chain, the two noninteracting spins attain effec- 
tive interaction between them. To study the dynamic 
detail we consider the system depicted in FiglTJ The 
Hamiltonian of this model reads 



where 



(1) 



(2) 



is the unperturbed Hamiltonian for the external spins, 

TV 

is the Hamiltonian of the transverse field Ising model, 
and 



N 



a=A,B j = l 



(4) 



is the interaction between the external spins and the en- 
vironment. Here cr^ and {a = A, B , f3 — x, y, z) are 
Pauli operators for the two external spins and the Ising 
chain respectively, cr^^- = (c^j + ic^ j)/2 are the corre- 
sponding raising and lowering operators, Ja/^fN is the 
homogeneous coupling constants between ath external 
spin and jth site of Ising chain with N being the number 
of sites in the chain. 

First of all. He is diagonalized with the combination 
of the Jordan- Wigner transformation (Tsj 



exp(7riE'^fef^fe)'^j 
fe=i 



(5) 



and the Bogliubov transformation 
1 ^ 

where for J/T^l, 

(7) 
(8) 



Okj = \l sm{kj) for A; > 0, 



okj = V ]y cos(fcj) for A; < 0, 



ipkj = - — [{I + Xcosk)(j)kj + Xsinkip^k]]- (9) 
Afc 



Here, 



Ak = Vl + A2 + 2Acosfc 



(10) 



is the energy spectrum of the quasi particle with k — 
2'Km/N, m = -N/2, • • • , N/2 - 1 for even iV, and m = 
~{N - l)/2, • • • , (A^ - l)/2 for odd N. For A = 1 and 
TO = -N/2, 

Afe = 0, (/)fe, = Vfcj = ±y^. (11) 

Thus, in the quasi-particle representation the Hamilto- 
nian of the TFIM is rewritten as 



HE = 2VY,kk[lfk^k-\). 



(12) 



with corresponding eigen state and energy being 

k 

Em = 2rEAfc"fe-rE^fc' (14) 

respectively. Here, |0) is the ground state and Uk = 'Hk'Hk 
is the particle number operator. 



III. EFFECTIVE HAMILTONIAN 

Generally speaking, Frohlich transformation [l^, [l^ is 
widely used in condense matter physics. It can solve a 
class of problems such as the induced effective interaction 
between two electrons by exchanging phonons with the 
crystal lattice. In this paper, the one-dimensional Ising 
chain plays the role as a medium to induce the effective 
interaction between the two external spins. Therefore, by 
virtue of Frohlich transformation, we obtain the effective 
Hamiltonian between the two external spins by tracing 
over the degrees of the environment. 

With an appropriate anti-Hermitian transformation S 
defined by the matrix elements 



(m| S \n) 



(to| Hi \n) 



(15) 
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which meets the condition Hj + [Hq,S] = 0, the effec- 
tive Hamihonian is approximated to the second order as 
Heff ^ He + Hei , where 



^ ^{0\[Hi,S]\0) 



(16) 



i ^((0| Hj \m) (m| S |0) - (0| S \m) (m| Hj |0)) 



with {|to)} and {Em} being the eigen states and eigen 
energies of He respectively, and |0) its ground state. 
The DSF [13 



N 



(17) 

(a, (3 = x,y) is the Fomier transformation of the Green 
function for the TFIM. It is calculated as 

„ „ (01 cr" I m) (ml erf, JO) 



n— 1 m 
N 



:5^e''="^(0k«|m) (mkf^jO) 



n—1 771 



x[p( 



Eq — E„i + Lu 



) - inSiEo - E„i + uj)] 



with p(l/x) being the principal value of 1/x. 

The right hand side of Eq. p6|l contains the following 
terms, which are expressed in terms of DSF's as 

5^ (0|cT;|m)(m|4|0)p(-— -1 



Eq — Em + UJ 

(0| oj \m) (m| CTj, |0) p(- ^^^) 



(19) 



E (Ok;|m)(mk,^|0)p( „ _ ^ ) 



N 

T 



-Bo — Em + 

{lm[5^^(0, w) + Syy{0, uj)] - 2Re5^^(0, uj)l?0) 



J2 (Ok7|m)(m|4|0)p( 



1 



Eq — Em + w 



^{lm[5^^(0, Lu) + Syy{0, Lu)] + 2ReS'^y(0, w)](21) 



Then, we rewrite the effective Hamiltonian as 



in terms of DSF's, where UM 



(22) 



MA = M + 4{Im[5""(0,M) + 5™(0,Ai)-5"^(0,-M) 

-syy{o, + 2Re[^-no, m) + s^y{o, -m)]}, 
f^B = ^i + Jl{lm[s^''{o,^l) + syy{o,^^)-s^^{o,-fI) 

_5W(0, -^)] + 2Re[^-^(0, m) + ^-''(0, 

91 = JAJB{Ms^''{o,^J-) + syy{o,^i) + s^^{o,-^l) 
+syyio, -/i)] + 2Re[5-^(o, m) - -/i)]}, 
52 = JAJBMS''''{o,^i)-syy{o,^i) + s^^o,-^l) 

^Syy{0,~^^)]. (23) 

In the forthcoming section, by using the numerical 
method in Ref.[l3], the DSF S'"^(fc,w) and thus the ma- 
trix elements of the effective Hamiltonian ([25)1 arc cal- 
culated explicitly. For further details about the diago- 
nalization method of the Ising-like model and the fast 
scheme for the calculation of Pfaffian, please refer to 

Refs. [n, E H nu, m m . 



IV. CRITICAL COUPLING 

In order to calculate the DSF numerically, we shall 
summarize the numerical method introduced in Ref . [l7| . 
For a spin chain in a transverse field with open ends, the 
Hamiltonian is described as [2J| 



H' 



N 



N-l 

E 



(24) 



After the Jordan- Wigner transformation, the Hamilto- 
nian is transformed into Fermion representation. Then, 
it is equivalent to solving the following eigen problem 

nam. 



$fc(A-B)(A + S) 



(25) 
(26) 



where A and B are two N x N matrixes with their 
matrix elements being Aij = 2fldij -t- JiS^+ij- + JSi-ij 
and Bij = JSi+ij — JSi^ij. According to the Wick- 
Bloch-de Dominicis theorem, the x-x correlation func- 
tion can be expressed in the form of the PfafSan of the 
2(2j + n — 1) X 2{2j + n — 1) antisymmetric matrix con- 
structed from elementary contractions 
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Pf 











j+n) 



j+nl 



-{ip+{t)ip+) -(ip^(t)ip+) -{Vt(t)^t) 



(27) 



where 



4 



4 ± 



(28) 



is the hnear comibination of the Fermion operators, and 

N 



N 

N 

p=l 
N 



(29) 
(30) 
(31) 
(32) 



p=i 



are the elementary contractions of zero temperature ob- 
tained from the finite temperature counterparts in Refs. 

The Pfafhan is the square root of the determinant of 
the corresponding antisymmetric matrix. A fast compu- 
tation scheme is given in Ref . [l9l. [23j . For an x iV 
antisymmetric matrix 



X 



A B 

-B^ C 



(33) 



with the dimensions of A, B, C being 2 x 2, 2 x (A^ — 2), 
(A^ — 2) X (A^ — 2) respectively, the Pfafhan of X can be 
computed in the following way. Because of 



h 

l-N-2 



X 


'h 


-A-^B 




'A 





In-2 




C + B'^A-^B 



with /„ being n dimensional unit matrix, wc have 

Det(A:) = Det(A)Det(C + B'^A-^B). (34) 
Since antisymmetric A is of the simple form 



A 



Xi2 
-Xi2 



(35) 



the N — 2 dimensional matrix C + B^A^^B is also anti- 
symmetric. The above procedure can be repeated times 
and times again. And the original matrix X is decom- 
posed into A^/2 2D antisymmetric matrices. Finally, due 




FIG. 2: (color online) Concurrence evolution vs A for A = 40, 
fj, = 1, Ja ~ Jb = 0.1 with an initial state \eg). 



to Pf(A) — xi2, the Pfafhan of matrix X will be simply 
a product of A^/2 numbers obtained from those 2x2 ma- 
trices in the above procedure. As it is not necessary for 
A to be a 2D matrix at the upper left corner X, A can 
be chosen to be a diagonal block such that Det(yl) is the 
largest, for the stability of the algorithm. 

Furthermore, other DSF's can be calculated according 
to the relation between correlation functions, that is 



(2f7)2 dt^ 



(36) 
(37) 



In the last section, we have obtained a typical spin-spin 
coupling in the effective Hamiltonian induced by the Ising 
chain. Driven by this Hamiltonian, two external spins 
can be entangled dynamically. To characterize the extent 
of entanglement, we use concurrence to measure the in- 
duced entanglement. For an arbitrary state of two-qubit 
system described by the density operator p, a measure of 
entanglement can be defined as the concurrence [1^, [2^ , 



C{p) = max{0, Ai — A2 — A3 — A4}, 



(38) 



where the A^ 's are the square roots of the eigenvalues of 
the non-Hermitian matrix pp in decreasing order. And 



(39) 
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FIG. 3: (color online) gi vs A for different A'^'s with /i = 1 
and Ja = Jb = 0.1. Red solid line for A'^ = 20, Red solid line 
for N = 20, green dashed line for A'^ = 30, blue dash-dot line 
hne for A'' = 40, black dotted hne line for A'^ = 50. 



FIG. 5: (color online) Fidelity vs A for different A'^'s with 
fi — 1 and Ja = Jb = 0.1. Blue solid line for A'^ = 4, black 
dashed line for A'' = 6, red dotted line for N = 8. 
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FIG. 4: (color online) g2 vs A for different A^'s with fi = 1 
and Ja = Jb = 0.1. Red solid line for A'' = 20, Red solid line 
for A'^ = 20, green dashed line for A'^ = 30, blue dash-dot line 
hne for A'' = 40, black dotted hne line for A^ = 50. 



where p* is the complex conjugate of p. 

As shown in Figl^l we investigate the evolution of 
concurrence under Heff. The two external spins start 
with an initial product state \eg). As time passes by, 
the two spins will eventually evolve into a maximum en- 
tangled state. It can be seen from the figure that the 
time needed for reaching maximum entanglement can be 
greatly shortened in the vicinity ofA = ri/J=l. In 
other words, the induced entanglement between two ex- 
ternal spins can be enhanced by quantum criticality. 

Furthermore, the relation between gi and A is plotted 
for different iV's in FigO Although the curves oscillates 
shapely, they share a characteristics in common. There is 
a peak near the critical point A = 1 for different curves. It 
is a reasonable result since it has been discovered that the 
entanglement between two nearest neighbors of the Ising 
chain achieves maximum near the critical point [l|, 
Similar result is also obtained for the relation between g2 



and A in Figdl However, besides the one near A = 1 there 
is another peak around A = 1.5 except that for iV = 40. 

We notice that in the above calculations, the thermo- 
dynamical limit condition is not used. As A varies from 
to 2, the single particle energy spectrum goes to con- 
tinuum as N increases. Since Frohlich transformation 
is equivalent to the second order perturbation theory, it 
may not be valid to apply it to obtain effective inter- 
action between the two spins. Therefore, we resort to 
mixed-state fidelity to demonstrate our approximation 
used in this paper. The mixed-state fidelity is given as 



F{pa,pi) = tr\/pl/^popy^ 



(40) 



which measures the degree of distinguishability between 
the two quantum states pq and pi. It has already been 
applied to the research on QPT [29]. Starting from the 
original Hamiltonian ([1]), we obtain the reduced density 
matrix po for the ground state of the two spins by first 
diagonalizing H and then tracing over all the degrees of 
the Ising chain. On the other hand, we can also obtain 
the density matrix pi of the ground state from the ef- 
fective Hamiltonian (|22|) of the two spins. In FiglH the 
relation between Fidelity and A is plotted for different 
TV's. Despite oscillations in some paramter intervals, the 
numerical method shows high fidelity over the whole in- 
terval. Moreover, curves of different TV's converge as A 
goes larger. 



V. CONCLUSION AND REMARK 

In summary, we have studied the dynamical process 
of two external spins simultaneously coupling to a trans- 
verse field Ising chain. With Frohlich transformation, we 
have deduced the effective Hamiltonian between these 
two spins. The matrix elements of the Hamiltonian are 
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expressed in terms of the DSF, which can be numeri- 
cally calculated. Through the numerical simulation, it 
is shown that the induced entanglement is enhanced by 
QPT. And the effective coupling constants reach maxi- 
mum near the critical point. By virtue of mixed-state 
fidelity, we demonstrate the validity of Frohlich transfor- 
mation and numerical simulation. Thus, the measure of 
the two spin entanglement can be an illustration of QPT. 

Besides the enhancement of the coupling intensity 
around the critical point, there are oscillations elsewhere. 
We remark that QPT takes place in the thermodynami- 
cal limit, i.e., N oo. In this limit, the energy spectrum 
of TFIM goes to continuum from zero to infinity at some 
parameter, i.e., A = 1. Therefore, the eigen energy of the 
external spins will definitely be resonant with one of the 
eigen energies of the Ising chain. In this circumstance, 
our method may not work well. We notice that a numer- 
ical method, the time-dependent density matrix rcnor- 



malization group (i-DMRG) [s^, was applied to central 
spin models, which were quite similar to ours. There- 
fore, in the near future, we will apply this method to 
our model to obtain a better result. We expect smooth 
curves similar to Fig. 1 in Ref. [l| . 
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